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Abstract 

I describe the first steps in the construction of semiclassical states for non-perturbative 
canonical quantum gravity using ideas from classical, Riemannian statistical geometry 
and results from quantum geometry of spin network states. In particular, I concentrate 
on how those techniques are applied to the construction of random spin networks, and 
the calculation of their contribution to areas and volumes. 

Introduction 

The non-perturbative canonical approach to quantum gravity has made a great deal of 
progress over the past decade, with a few of the main developments being a rigorous con- 
struction of the kinematical Hilbert space in the connection representation, results on the 
quantum geometry of states in this Hilbert space, and the proposal of a Hamiltonian oper- 
ator for the theory Q. Although the quantum theory is not yet complete even in the pure 
gravity case, enough of it is in place for us to start developing rigorously the semiclassical 
theory and its coupling to matter, with the aim of addressing its phenomenological aspects, 
some of which have already been treated in the literature 0, but only in a heuristic way. 

A basic ingredient for such developments is an understanding of semiclassical gravita- 
tional states and their geometry. At the heuristic level, these considerations led several 
years ago to the proposal of weave states 0. Starting from loop states ||, which encode 
distributional information about the spatial geometry along the one- dimensional submani- 
folds defined by a set of (possibly intersecting) loops, weave states were based on uniformly 
distributed loops in the spatial manifold, and sought to encode in some sense all of the 
geometry down to the length scale associated to the loop spacing. 

Analogs of the weave states in the context of the rigorous formulation of the quantum the- 
ory are natural candidates for semiclassical states, and some papers have recently appeared 
looking at this question ||. Among the first questions about such states that need to be 
addressed are their definition and relation to the kinematical Hilbert space, and the nature 
and extent of the geometrical information they encode. These issues are being studied in a 
collaboration with A. Ashtekar ||, and I will describe here some of the work, touching on the 
more geometrical aspects of our proposal, and using as the main example a flat geometry. 

The basic, canonically conjugate variables of the theory are a density- weighted triad field 

on the spatial manifold S (the "electric field" ) , from which the inverse of the spatial metric 
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q a b can be obtained as (det q) q ab = EfE h - Sy, and an su(2)-valued connection A % a = Y' l a + kK % a1 
where T l a is the connection with respect to which Ef is covariantly constant, K l a corresponds 
to the extrinsic curvature K a b of a constant t slice on a classical solution, in the sense that 
K\ = (det q)' 1/2 K ab E\, and k is an arbitrary number, the Immirzi parameter. 

In the connection representation, pure states in the kinematical Hilbert space 7i° of 
the theory (i.e., before imposing the quantum constraints) are given by functionals ^f[A] of 
suitably generalized connections A on E, and a basis for Ti° is given by the spin network 
states 0, defined by triples S = {7, j, I}, where 7 is a directed graph embedded in E, j = {j e } 
a set of labels for representations of SU(2), one for each edge e of 7, and I = {I v } a set of 
labels for intertwiners at all vertices v, i.e., a (gauge-invariant) prescription for contracting 
the indices of the matrices RP e ^ that the representations j e associate with the holonomies 
U(e, A) of A along the edges e incident at each v. For every such triple, a wave functional 
in the connection representation for quantum gravity is defined by 

= (1,3,1 I A) := J{Jl v I v R^\U{e,A)) . (1) 

On these states, of course, A % a is a multiplication operator, Ef a differential one. 

In terms of (p]), weave states have a graph 7 given by a uniformly distributed collection of 
loops, with respect to a classical geometry (E, q a b), and all holonomies are taken in the j = \ 
representation. For these states, it was found that they were eigenstates of the operators 
A s and Vr corresponding to areas of surfaces S and volumes of regions R in E, as is true 
more generally for spin network states, and that the eigenvalues would be the correct areas 
and volumes given by the classical metric, for large surfaces and regions, if the loops were 
uniformly scattered with mean spacing a = ^/2n £p between them, with £p the Planck length, 
which was seen as a first indication of a fundamental discreteness in the non-perturbative 
quantum theory, emerging from the theory itself as opposed to being put in by hand. 

In general, in order to ask whether a quantum state \l/ is semiclassical, a set of observables 
{g a } must have been previously specified; the state is semiclassical if the expectation values 
of the corresponding set of operators coincide with the values assigned by those observables 
to a classical phase space point (E,A) on E, (g a )^ = g Q (E,A), and the uncertainties are 
small, in the sense that (Ag a )^ <C g a (E,A), or (g 2 )^ <C \g a {E, A)] 2 (except when the 
classical value itself is small). A proposal for a set of semiclassical states will then start with 
the choice of a sufficiently large, and physically relevant, set {g a }. 

In the weave state approach, the choice was {g a } = {A s , Vr} for large and "slowly 
varying" S and R, and we will make the same choice here. Other choices are possible, 
and one is illustrated by the coherent state construction proposed by T. Thiemann and 
collaborators ||, in which the basic observables are the holonomies of A along the edges of 
a given graph 7, and what can be viewed as fluxes of the electric field through a given set 
of surfaces, each one intersecting transversally one of the edges of the graph. 

In this paper, we give a procedure for constructing random, uniformly distributed graphs 
in E that can be used to define random spin networks. Although the procedure is more 
generally applicable, we treat the case of a flat metric, and use results from classical statistical 
geometry of Euclidean manifolds and from quantum geometry of spin network states to find 
the contribution to areas and volumes from states of this type; such states are eigenstates of 
those geometrical operators, and can be thought of as corresponding to the previous, heuristic 
weave states. We then comment on the possibility of using those states to construct more 
general ones, and on the extension to curved geometries. 
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Classical Statistical Geometry: Random Complexes 

Consider a Euclidean manifold (E, e a j,), of finite volume V% and arbitrary dimension D, for 
the time being. We start by sprinkling N points at random in (£,e a b), independently and 
with uniform density p = iV/V^; that is, the probability density that each sprinkled point 
fall at any location x in £ is 



dv \ /em 

dP(x) = — = ^±ld D x, (2) 

and the probability that it fall in any (measurable) region R C S is therefore Vr/Vs (all 
volumes Vr = V(R) will be defined using dv = yfed D x from now on.) Combining these prob- 
abilities for single points we find that, when the N points have been chosen, the probability 
that exactly n of them be in any given R C S is given by the distribution 

PbinomiaKn, J2|J\T, S) = ( N ) (^] ( 1 - Y^)" " . (3) 



To perform an actual sprinkling in a computer simuation is very easy if the geometry is flat, 
since then coordinate values can be chosen uniformly at random in a Cartesian chart; if the 
geometry is curved, other well-known techniques can be applied (see, e.g., Ref. 0). 

Two special limiting cases deserve mention. One is the continuum limit, approached as 
iV and p become very large, with constant Vs; the other is the infinite volume limit, in which 
N and Vs are very large, with p = N/Vy, a constant. Let us analyze the latter situation in 
more detail. Since V% = oo, in this case we cannot use the probability density (Q); we do, 
however, have probabilities for finite regions. It is a standard, well-known result, that in 
the infinite volume limit the distribution (|3J) approaches a Poisson distribution with mean 
n = p V R , 

Pbinomial(n, R\N, S) « P PoissoI1 (fl | P V R ) = J (4) 

for this reason, uniform random distributions are often called Poisson random lattices [ [Tcfl . 
Thus, simulating a random sprinkling in a region R which is part of an infinite-volume 
manifold S is a two-step process, in which one first generates the number Nr of points 
in R using the distribution -Ppoisson^Rlp^i?,), an d then generates locations for those points 
inside R as in the compact manifold case. Of course, the same two-step procedure can be 
followed when a compact manifold needs to be divided into two or more disjoint regions for 
the simulation, although in that case the binomial distribution is used for the first step. 

Regarding fluctuations in n, the standard deviation of a Poisson distribution is well- 
known, <jP° lsson = = yj pVr] for the binomial distribution, one can readily verify that 



^binomial / „t / ( i P^R \ ^.Poisson ( i P^R \ / r \ 

where in the last, large N approximation we have neglected terms of order iV -2 . Similar 
results hold for individual probabilities, which justifies the use of the Poisson distribution 
as an approximation even when the number iV of points in a compact manifold is fixed; one 
may say, e.g., that the probability that the region R contain no sprinkled points is e~ pVa , 
using Eq. (§), provided N is large. 
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Once an arbitrary, locally finite (assumed here to be random) distribution of points {p,} 
is given in a Euclidean manifold (E, e Q &), there is a well-known construction, often called 
Dirichlet-Voronoi construction, of two cell complexes based on that distribution: a simplicial 
complex A triangulating E, with those points for vertices, and its dual cell complex Q. 

The simplicial complex A can be obtained in the following way. We associate a simplex 
with a subset of D + 1 points among all the {p^} if the (unique) S D_1 sphere passing through 
all of them does not contain any other point pj (this sphere is defined by the metric e a b, 
and has as its center the unique point which is equidistant from the D + 1 p^s). The set of 
all such simplices is the desired complex; it covers E, and any two of them can only have a 



vertex, an edge, or a 2-face in common |TT]]. Being made of simplices, this complex has a 
fixed number of /-faces in each fc-face, for / < k (for example, in D = 3 dimensions, there 
are exactly 4 triangles, / = 2, in each tetrahedron, k = 3), but the number of /-faces sharing 
a given fc-face, for / > k, depends in most cases on the specific set of points used, and for us 
will be a random variable (like, e.g., again for D = 3, the number of edges, / = 1, sharing a 
given vertex, k = 0). 

The cell complex fl, dual to A, is obtained by defining, for each point p iy a cell to be 
the set of all manifold points which are closer to pi than to any other pj; the set of all such 



cells as pi varies is the desired complex |Q|. Therefore, each D-cel\ Qi of this complex is 
dual to a vertex pi of the simplicial complex, and each (D — l)-cell fly is dual to an edge pij, 
and perpendicular to it if they meet.Q In general, there is a duality between /c-faces Pi x ...i h in 
the simplicial complex and (D — k + l)-cells ifc , and the incidence relations reflect this 
duality. Thus, there is a fixed number of /-cells sharing each fc-cell, for / > k (for example, 
in D = 3 dimensions, there are exactly 4 edges sharing each vertex, except for degenerate 
cases), but the number of /-cells in each k-ce\\, for / < k, in most cases depends on the 
specific set of points used, and for us will be a random variable (like, e.g., again for D = 3, 
the number of faces in each 3-cell). We denote the set of /c-faces of A by A^ k \ and the set 
of /-cells in fl by fl®. 

When the above constructions are used with uniform random distributions of points of 
density p, one gets random simplicial and dual cell complexes A p and fl p , that have been 
studied for a long time. In the context of gauge theory, their use was proposed in the early 
80's as a way of implementing a short distance cutoff without breaking Euclidean invari- 
ance [12|, [11], but results on statistical properties of random complexes had been obtained 
long before in metallurgy and mineralogy fjH|, motivated by studies of crystal formation by 
random nucleation in minerals, and by mathematicians |L4] . 

For our applications to quantum gravitational states, two kinds of related properties of 
the above complexes will be important. Incidence relations between simplices or cells of 
different dimensionalities, that are statistical topological properties of the complexes, will 
relate the number of D-simplices to the number of sprinkled points, and therefore give us 
the density of dual cell complex vertices in terms of p. On the other hand, metric properties, 
such as average cell sizes, will be more closely related to average intersection numbers of cells 

1 With the above, standard construction of A and fi, not all edges meet their dual faces. One possible 
slight modification of the construction is to use, given a set of sprinkled points, the same simplicial complex 
A, but choose as vertices of the dual complex the incenters (centers of the inscribed spheres) of all simplices, 
as opposed to the circumcenters (centers of the circumscribed spheres) , which is what the standard procedure 
amounts to. The topology of the resulting dual complex £1' is the same as that of f2, but in the pair (A, f2') 
dual elements always intersect, although in general not perpendicularly. 
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of various dimensionalities with given subsets of S, such as the number of intersections per 
unit area between a given surface and graph edges. All properties of this type are known 
for two- and three-dimensional Euclidean space; we will discuss here those properties we will 
need later for our quantum statistical geometry results. 

To calculate incidence relations, we consider the (finite) simplicial complex A p . Denote 
by Nk the total number of fc-faces of A p ; in particular, N = A" is the number of points. Two 
relations among the AVs can be used in any dimension. One is the expression for the Euler 
characteristic of a simplicial complex, Y^k=o(~ ^) k Nk = x(A p ), where the dependence on the 
global topology of S can be eliminated by dividing the equation by N, and then taking the 
large A" limit, in which x(A p )/A^ — > 0, and 

The other relation is the fact that each D-simplex has D + 1 faces of codimension 1, each 
one shared by two D-simplexes, so that 

iV.-i = ^ N D . (7) 

In D = 2 dimensions, these two equations would be sufficient for determining the average 
number of edges and 2-simplices per unit volume, which in the uniform sprinkling case are the 
mean densities of the corresponding simplices, constant throughout E; in D = 3 dimensions, 
however, we need more equations. A detailed calculation, where one integrates explicitly 
over the probabilities of finding points at various locations (see, e.g., Ref. fllj), shows that 
each cell has on average 967r 2 /35 vertices; since each of those vertices is shared by four cells, 
and there are p cells per unit volume, we find the density of cell complex vertices to be 

N 3 1 96tt 2 24tt 2 

% = 1^ P= ^ P - (8) 

To calculate the metric properties of random complexes, we must make extensive use of 
the probabilities associated to the point distribution. In D = 3 dimensions, we are interested 
in the mean number of edges intersected by a flat surface per unit area, for which we can 
use the following argument. As has long been known (see, e.g., references in []13||), from 
a generalization of the Buffon needle method of calculating it by random tosses of a stick 
on a series of parallel lines, when a randomly oriented set of lines of arbitrary (possibly 
disconnected) shape and total length L, contained in a region of volume V, is cut by a 
surface of area A, the number of intersections between the lines and the surface per unit 
area is very simply related to the line length per unit volume, N int /A = | L/V . The problem 
of finding the mean number of intersected edges is thus reduced to that of finding the mean 
total edge length per unit volume. A known calculation gives that the mean edge length per 
3-cell is (47r) 5//3 (3 1 / 3 /5) T(|) p -1 ^ 3 , from which L/V can be found multiplying by the mean 
number of cells per unit volume, p, and dividing by the number of cells sharing each edge, 
3. Putting all of this together then gives 

(4tt)^ 3 _3 1 / 3 
30 

This result does not depend on the shape of the surface (in particular, S does not have to 
be flat), but other moments of the probability distribution for N int , such as its width, do. 



N int /A = v J on T(|) p 2 / 3 « 2.917p 2 / 3 . (9) 
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Quantum Statistical Geometry: Areas and Volumes 

In this section, we start discussing how to use the random cell complexes introduced above 
to construct quantum gravity states. The amount of structure we use from each pair of 
complexes is a choice we make in each approach, depending largely on the set {g a } of 
observables we want to reproduce. Since our observables here are areas and volumes, results 
from quantum geometry that we will summarize below allow us to make the simplest possible 
choice: We use only the 1-skeleton, or set 7 X = Cl^ of edges of the dual cell complex obtained 
with a randomly generated set of points x = {xj}, as a random spin network graph. The 
main advantage of this choice, as opposed to using the simplicial complex A p , is the fact 
that almost every vertex is exactly (D + l)-valent in D dimensions; notice that, given that 
the edges are geodesic segments, the degree of differentiability of the graph depends on the 
geometry, and in the present case the graph is piecewise analytic. 

In this paper, we will also make other simplifying choices, namely those of using only one 
such graph at a time for the states we consider, and of defining a spin network state ^ 7 j,i 
with the same spin label j at each edge and the same intertwiner I at each vertex; we can 
then write our states as ^ 7 ,j,/. These further choices are made purely for illustrative purposes 
here, and will be changed in a more complete treatment; a few comments on this point will 
be included in the final section, and see Ref. || for details. This amounts to assuming that 
any single random graph 7 provides a good sampling of the underlying manifold, and we end 
up with a two-parameter family of states for each such graph. One question we can address 
is then what constraints are placed on the parameters by the requirement that the area and 
volume operators have the right values on these states. 

Consider a state \^ 7 j,j of the type just described. Given a surface S in S, any given 
vertex of the random graph 7 will fall on S with probability zero, and we can consider all 
intersections of S and 7 to be single, transversal edges. In this case the spin network is an 
eigenstate of the operator with eigenvalue [TH] 

A S =87Tk£lJ2 a \/Uja + l), (10) 

where j a is the half-integer label for the a-th edge crossing 5*. If all the j a 's are equal to a 



given j, the area eigenvalue for a given spin network becomes Sn k £p N int ^j j(j + 1), where 
iVj nt is the number of intersections between S and the graph, so from Eq. @ we obtain that 
on average classical and quantum areas agree if p and j satisfy 

( 47r )5/3 3 l/3 
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T(l)p^=(87rkJjU^T)eiy . (11) 



Given a region R C £, spin network states are eigenvectors of the volume operator Vr, 
and the corresponding eigenvalues receive a contribution from each vertex of the graph, 

V R % jjtI = k (8vrfc) 3/2 4 J2 V \f\&\%,j,l , (12) 

where the constant k,q is an undetermined factor arising from a regularization ambiguity for 
the volume operator, and q v is an operator corresponding to the determinant of the spatial 
metric at v , whose eigenvalues are determined by the j a 's of all edges incident at that vertex 
and the intertwiner I v While a closed formula like the one for area eigenvalues, Eq. 



(pi]), is not available for volumes, the calculations simplify in the case of four-valent vertices. 
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When the four j's are equal, the number of independent intertwiners at v, and thus of 
eigenvalues of q v , is j + 1; when j is even, one of those eigenvalues is zero, while the others, 
and all eigenvalues when j is odd, come in pairs of opposite signs and therefore give doubly 
degenerate volume eigenvalues. The volume eigenvalues have been calculated for the first 
few values of j |T7|]; for example, the only eigenvalue for j = 1, and the non-zero ones for 



j = 2, are 

ol/4 ol/4 

Ai = ^ r «ofc 8/2 4, A 2 = ^^ 3/2 4- (13) 
In general, we can write the eigenvalues for vertices of our type in the form 

\ jt i = f(j,I)Kok* /2 4, (14) 
so from Eq. (H) we obtain that on average classical and quantum volumes agree if 



24tt 2 
~35 



P 



(f(jJ)K Q 0Hl)~ l . (15) 



If the parameters k and kq of the theory are known, equations (|Tl|) and (|Tj|) are constraints 
on the parameters characterizing the random spin network states. Notice, however, that 
more constraints on the latter parameters are expected to arise from matching the manifold 
curvature in more general cases, in which we consider non-flat manifolds (£, q a b)i and we can 
also think of those two equations as fixing k once values for j and / have been determined, 

(47r) 8/3 gl/3 X 3/2 247r 2 



j(j + i)^Y 5 — r(f) =_/(;,/)«„. (i6) 



Concluding Remarks 

An obvious way to generalize the states used in this paper is to relax the assumption that 
all j's and all Fs on a given graph 7 be equal, and define a state 

l7,C^=£ w C7 Jf i|7j,I>, (17) 

for some set of coefficients C- } \ to be determined. But the more interesting modification to 
our states comes from the observation that, given (£, e^), the sprinkling process does not 
give a unique set of points in S, but rather a probability density for iV-point distributions, 



dP p ( Xl , x N ) = M J] ^ d° Xl , (18) 

i=l y z 

parametrized by the density p, or N = pV^, and therefore the dual cell complex construction 
also gives us a probability density on the set of graphs embedded in S depending on the 
parameter p, rather than a single graph. We can then use this probability density to integrate 
over random graphs, and obtain either pure states resulting from their superposition, or 
mixed states of the form 

* e (p,C) = J dP p (x 1 ,...,x N )\^,C)(^,C\ , (19) 
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where the subscript e makes the dependence on the metric e a b explicit. Superposing graphs 
in this way may seem like a complication, but on the other hand it has the advantage that, 
while a single |7 x ,j,I) would only have an approximate Euclidean invariance, the integral is 
exactly invariant (despite having an ultraviolet cutoff scale), as well as covariant with respect 
to the action of diffeomorphisms. 

In the previous section, we have only discussed (for the simpler states ^-yjj) the expecta- 
tion value of areas and volumes, and seen how they give rise to conditions on the parameters 
the states depend on. The uncertainties that states of the form (|l^) associate with those 
observables, or others we may use to identify semiclassical states, will have a quantum con- 
tribution and a classical, statistical one; imposing that they be smaller than the desired 
tolerance will introduce further conditions, and restrictions on the length scales defined by 
the geometry and the subsets of £ we consider. 

In addition, it is useful to keep in mind alternatives to some other choices we made; 
although they were the simplest ones, there is no guarantee that once we understand the 
dynamical aspects better they will appear as the best ones. Specifically, in a curved manifold, 
uniform distributions of points are not the only possible covariantly defined ones; a uniform 
distribution in a Riemannian manifold (E, q a b) is equivalent to using the density dv = ^/qd D x 
as a measure dfi on E, while alternative ones can be defined using any scalar constructed 
from the metric, such as the Ricci scalar, dfi = Rdv, or any other curvature scalar. Possibly 
related to this is the fact that a better understanding of the theory may show that, in order 
for a spin network state to have a fully consistent semiclassical interpretation (including, 
e.g., the fact of giving rise to distributions for the spatial geometry and its time derivative 
which are both peaked around classical values) , the weave we use must be constructed using 
correlations in the sprinkling process, that were not used here. 
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